M. T. Schobeiri 
B. Oztiirk 

Turbomachinery Performance and 
Flow Research Laboratory (TPFL), 
Texas A&M University, 
College Station, TX 77843-3123 


David E. Ashpis 

National Aeronautics and Space Administration, 
John H. Glenn Research Center at Lewis Field, 
Cleveland, OH 44135 


Effect of Reynolds Number and 
Periodic Unsteady Wake Flow 
Condition on Boundary Layer 
Development, Separation, and 
Intermittency Behavior Along the 
Suction Surface of a Low 
Pressure Turbine Blade 

The paper experimentally studies the effects of periodic unsteady wake flow and Reynolds 
number on boundary layer development , separation , reattachment , and the intermittency 
behavior along the suction surface of a low pressure turbine blade. Extensive unsteady 
boundary layer experiments were carried out at Reynolds numbers of 110,000 and 
150,000 based on suction surface length and exit velocity. One steady and two different 
unsteady inlet flow conditions with the corresponding passing frequencies, wake veloci- 
ties, and turbulence intensities were investigated. The analysis of the experimental data 
reveals details of boundary layer separation dynamics which is essential for understand- 
ing the physics of the separation phenomenon under periodic unsteady wake flow and 
different Reynolds numbers. To provide a complete picture of the transition process and 
separation dynamics, extensive intermittency analysis was conducted. Ensemble- 
averaged maximum and minimum intermittency functions were determined, leading to the 
relative intermittency function. In addition, the detailed intermittency analysis was aimed 
at answering the question as to whether the relative intermittency of a separated flow 
fulfills the universality criterion. [DOI: 10.1115/1.2219762] 


Introduction 

In recent years, gas turbine engine aerodynamicists have fo- 
cused their attention on improving the efficiency and performance 
of the low pressure turbine (LPT) component. Research per- 
formed by the industry, research centers, and academia has shown 
that a reduction of the blade number can be achieved without 
substantially sacrificing the efficiency of the LPT blading. This 
reduction contributes to an increase in thrust/weight ratio, thus 
reducing the fuel consumption. Unlike the high pressure turbine 
that operates in a relatively high Reynolds number environment, 
the LPT of large commercial engines operates at Reynolds num- 
bers ranging from 75,000 to 400,000. Since the major portion of 
the boundary layer, particularly along the suction surface, is lami- 
nar, the low Reynolds number, in conjunction with the local ad- 
verse pressure gradient, makes it susceptible to flow separation, 
thus increasing the complexity of the LPT boundary layer aerody- 
namics. The periodic unsteady nature of the incoming flow asso- 
ciated with wakes that originate from upstream blades substan- 
tially influences the boundary layer development, including the 
onset of the laminar separation, the extent of the separation bubble 
and its turbulent reattachment. Of particular relevance in the con- 
text of LPT aerodynamics is the interaction of the wake flow with 
the suction surface separation bubble. While the phenomenon of 
the unsteady boundary layer development and transition in the 
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absence of the separation bubbles has been the subject of intensive 
research, the multiple effects of mutually interacting parameters 
on the LPT boundary layer separation and their physics still re- 
quires more research for full understanding. 

The significance of the unsteady flow effect on efficiency and 
performance of compressor and turbine stages was recognized in 
the early seventies by several researchers. Fundamental research 
by Pfeil and Herbst [1], Pfeil et al. [2], and Orth [3], studied and 
quantified the effect of unsteady wake flow on the boundary layer 
transition along flat plates. Schobeiri and his coworkers [4-7] 
experimentally investigated the effects of the periodic unsteady 
wake flow and pressure gradient on boundary layer transition and 
heat transfer along the concave surface of a constant curvature 
plate. The measurements were systematically performed under 
different pressure gradients and unsteady wake frequencies using 
a squirrel cage type wake generator positioned upstream of the 
curved plate. Liu and Rodi [8] carried out the boundary layer and 
heat transfer measurements on a turbine cascade, which was in- 
stalled downstream of a squirrel cage type wake generator, men- 
tioned previously. 

Analyzing the velocity and the turbulence structure of the im- 
pinging wakes and their interaction with the boundary layer, 
Chakka and Schobeiri [7] developed an intermittency -based un- 
steady boundary layer transition model. The analysis revealed a 
universal pattern for the relative intermittency function for all the 
frequencies and pressure gradients investigated. However, the 
above investigations were not sufficient to draw any conclusion 
with regard to an eventual universal character of the relative in- 
termittency function. Further detailed investigations of the un- 
steady boundary layer on a high Reynolds number turbine cascade 
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by Schobeiri et al. [9,10] and its subsequent analysis [11,12] veri- 
fied the universal character of the relative intermittency function. 
For this purpose, Schobeiri et al. [9] utilized a conceptually dif- 
ferent type wake generator, which is also used for the investiga- 
tion presented in this paper. Fottner and his coworkers [13,14] and 
Schulte and Hodson [15] used the same wake generating concept 
for the investigations on the influence of the unsteady wake flow 
on the LPT-boundary layer. Kaszeta, Simon, and Ashpis [16] ex- 
perimentally investigated the laminar-turbulent transition within a 
channel with the two curved walls resembling the suction and 
pressure surfaces of an LPT-blade using a retractable wake gen- 
erator. Lou and Hourmouziadis [17] investigated the effect of os- 
cillating inlet flow conditions on laminar boundary layer separa- 
tion along a flat plate under LPT-pressure conditions. This was 
emulated by contouring the test section top wall. They studied the 
Reynolds number effect on the transition region. Their results 
showed that the higher Reynolds numbers cause an earlier transi- 
tion and reduction of the transition length, while the separation 
point does not change its location. 

Using the surface-mounted hot film measurement technique, 
Fottner and his coworkers [13,14], Schroder [18], and Haueisen, 
Hennecke, and Schroder [19] documented strong interaction be- 
tween the wakes and the suction surface separation bubble on the 
LPT blades, both in the wind tunnel cascade tests and in a turbine 
rig. Furthermore, they investigated the boundary layer transition 
under the influence of the periodic wakes along the LPT surface 
and found that the interaction of the wake with the boundary layer 
greatly affects the loss generation. The investigations by Halstead 
et al. [20] on a large scale LP turbine use surface-mounted hot 
films to acquire detailed information about the quasi-shear stress 
directly on the blade surface. Investigations by Cardamone et al. 
[14] and Schroder [18] indicate that the benefit of the wake- 
boundary layer interaction can be used for the design procedure of 
modem gas turbine engines with a reduced LPT blade number, 
without altering the stage efficiency. 

Most of the studies mentioned above on LP turbine cascade 
aerodynamics have largely concentrated on the measurement of 
the signals stemming from hot films mounted on the suction and 
pressure surfaces of the blades under investigation. Although this 
technique is qualitatively reflecting the interaction of the unsteady 
wake with the boundary layer, because of the lack of an appropri- 
ate calibration method, it is not capable of quantifying the surface 
properties such as the wall shear stress. The few boundary layer 
measurements are not comprehensive enough to provide any con- 
clusive evidence for interpretation of the boundary layer transition 
and separation processes and their direct impact on profile loss, 
which is a critical parameter for blade design. Furthermore, the 
numerical simulation of the unsteady LPT blade aerodynamics, 
using conventional turbulence and transition models, fails if it is 
applied to low Reynolds number cases. Recent work presented by 
Cardamone et al. [14] shows that in the steady state case at Re 
=60,000 the separation is captured, however, for the unsteady 
case the separation bubble is not reproduced. 

A recent experimental study by Schobeiri and Oztiirk [21,22] 
investigated the physics of the inception, onset and extent of the 
separation bubble along a low pressure turbine blade, which was 
the first part of a series of investigations carried out at TPFL. A 
detailed experimental study on the behavior of the separation 
bubble on the suction surface of a highly loaded LPT blade under 
a periodic unsteady wake flow was presented in Ref. [21]. Surface 
pressure measurements were performed at Re =50,000, 75,000, 
100,000, and 125,000. Increasing the Reynolds number has re- 
sulted in no major changes to the surface pressure distribution. 
They concluded that the unsteady wake flow, with its highly tur- 
bulent vortical core over the separation region, caused a periodic 
contraction and expansion of the separation bubble. It was pro- 
posed that, in conjunction with the pressure gradient and periodic 
wakes, the temporal gradient of the turbulence fluctuation, or 
more precisely the fluctuation acceleration dv^J dt, provides a 
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higher momentum and energy transfer into the boundary layer, 
energizing the separation bubble and causing it to partially or 
entirely disappear. They found that for dv rms ldt> 0 the separation 
bubble starts to contract, whereas for £?r rms /f%<0 it gradually 
assumes the shape before the contraction. They argued that not 
only the existence of higher turbulence fluctuations expressed in 
terms of higher turbulence intensity influences the flow separa- 
tion, but also its gradient is of crucial importance in suppressing 
or preventing the onset and the extent of the separation bubble. 
They stated that the fluctuation gradient is an inherent feature of 
the incoming periodic wake flow, and does not exist in a statisti- 
cally steady flow that might have a high turbulence intensity. They 
also stated that unsteady wake flow, with its highly turbulent vor- 
tical core passing over the separation region, caused a periodic 
contraction and expansion of the separation bubble and a reduc- 
tion of the separation bubble height. Increasing the passing fre- 
quency associated with a higher turbulence intensity further re- 
duced the separation bubble height [21]. 

The objective of the present study, dealing with the specific 
issues of the LPT boundary layer aerodynamics, is to provide 
detailed unsteady boundary layer flow information to understand 
the underlying physics of the inception, onset, and extension of 
the separation bubble for different Reynolds numbers. Further- 
more, the unsteady boundary layer data from the present and 
planned experimental investigations will serve to extend the inter- 
mittency unsteady boundary layer transition model developed by 
Schobeiri and his coworkers [7,11,12] to the LPT cases, where a 
separation occurs on the suction surface at a low Reynolds num- 
ber at the design and off-design incidence. The experimental re- 
sults are also intended to serve as benchmark data for a compari- 
son with numerical computation using dns or rans codes. 

It is well known that the boundary layer measurement is one of 
the most time consuming aerodynamic measurements. Any at- 
tempt to increase the number of parameters to be studied would 
inevitably result in a substantial increase of the measurement time. 
Considering this fact, the research facility described in Refs. 
[9,10], with state-of-the-art instrumentation, has been substan- 
tially modified to systematically and efficiently study the influence 
of the periodic unsteady and highly turbulent flow on the LPT 
cascade aerodynamics at the design and off-design incidence 
angles, where the Reynolds number, wake impingement fre- 
quency, free-stream turbulence, and the blade solidity can be var- 
ied independently. 

Experimental Research Facility 

To investigate the effect of unsteady wake flow on turbine and 
compressor cascade aerodynamics, particularly on unsteady 
boundary layer transition, a multipurpose, large-scale subsonic re- 
search facility was designed and has been taken into operation 
since 1993. Since the facility in its original configuration is de- 
scribed in Refs. [9,10,21], only a brief description of the modifi- 
cations and the main components is given below. The research 
facility consists of a large centrifugal compressor, a diffuser, a 
settling chamber, a nozzle, an unsteady wake generator, and a 
turbine cascade test section as shown in Fig. 1. The compressor, 
with a volumetric flow rate of 15 m 3 /s, is capable of generating a 
maximum mean velocity of 100 m/s at the test section inlet. The 
settling chamber consists of five screens and one honeycomb flow 
straightener to control the uniformity of the flow. 

Two-dimensional periodic unsteady inlet flow is simulated by 
the translational motion of a wake generator (see Fig. 1), with a 
series of cylindrical rods attached to two parallel operating timing 
belts driven by an electric motor. To simulate the wake width and 
spacing that stem from the trailing edge of rotor blades, the diam- 
eter and number of rods can be varied. The rod diameter, its dis- 
tance from the LPT blade leading edge, the wake width and the 
corresponding drag coefficient are chosen according to the criteria 
outlined by Schobeiri et al. [23]. The belt-pulley system is driven 
by an electric motor and a frequency controller. The wake-passing 
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Fig. 1 Turbine cascade research facility with the components and the adjustable 
test section 


frequency is monitored by a fiber-optic sensor. The sensor also 
serves as the triggering mechanism for data transfer and its ini- 
tialization, which is required for ensemble- averaging. This type of 
wake generator produces clean, two-dimensional wakes, whose 
turbulence structure, decay and development are, to a great extent, 
predictable [23]. The unsteady boundary layer transition and heat 
transfer investigations [9-12] performed on this facility serve as 
the benchmark data for validation of turbulence models, transition 
models and general code assessments. 

To account for a high flow deflection of the LPT cascade, the 
entire wake generator and test section unit, including the travers- 
ing system, was modified to allow a precise angle adjustment of 
the cascade relative to the incoming flow. This is done by a hy- 
draulic platform, which simultaneously lifts and rotates the wake 
generator and test section unit. The unit is then attached to the 
tunnel exit nozzle with an angular accuracy less than 0.05 deg, 
which is measured electronically. 

The special design of the facility and the length of the belts 
(L belt =4960 mm) enable a considerable reduction of the measure- 
ment time. For the present investigation, two clusters of rods with 
constant diameters of d=2 mm are attached to the belts as shown 
in Fig. 2. The two clusters, with spacings 5^=160 mm and S R 
= 80 mm, are separated by a distance which does not have any 
rods, thus simulating steady state case (S R =™). Thus, it is possible 
to measure sequentially the effect of three different spacings at a 
single boundary layer point. To clearly define the influence do- 
main of each individual cluster with the other one, the clusters are 
arranged with a certain distance between each other. Using the 
triggering system mentioned above and a continuous data acqui- 
sition, the buffer zones between the data clusters are clearly 
visible. 

The data analysis program cuts the buffer zones and evaluates 
the data pertaining to each cluster. Comprehensive preliminary 
measurements were carried out to make sure that the data were 
exactly identical to those when the entire belt length was attached 
with rods of constant spacing, which corresponded to each indi- 
vidual cluster spacing. The cascade test section shown in Fig. 1, 
located downstream of the wake generator, includes five LPT 
blades with a height of 200.0 mm and the chord of 203.44 mm. 
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For boundary layer investigations, five identical “Pak B” airfoils 
designed by Pratt & Whitney were implemented and whose cas- 
cade geometry is given in Table 1 . 

The blade geometry resembles the essential feature, such as the 
laminar boundary layer separation, that is inherent to typical LPT 
blades. The blade geometry was made available to NASA re- 
searchers and academia to study the specific problems of LPT 
flow separation, its passive and active control and its prevention. 
As shown in Ref. [9], this blade number is necessary and suffi- 
cient to secure a spatial periodicity for the cascade flow. Identical 
pressure distributions taken from blade Nos. 2 and 4 (counted 
from the test section bottom, Fig. 1) demonstrate the required 
periodicity. These blades were specially manufactured for mea- 
surement of pressure and showed identical pressure distributions. 



t(m/s) 

Fig. 2 Wake generator 
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Table 1 Parameters of turbine cascade test section 


Parameters 

Values 

Parameters 

Values 

Inlet velocity 

V in =4 m/s 

Inlet turbulence intensity 

Tu in = 1.9% 

Rod translational speed 

U=5.0 m/s 

Blade Re number 

Re= 110,000 

Nozzle width 

W= 200.0 mm 

Blade height 

h B = 200 mm 

Blade chord 

c= 203.44 mm 

Cascade solidity 

a= 1.248 

Blade axial chord 

c ax = 182.85 mm 

Zweifel coefficient 

in 

(N 

II 

* 

Blade suction surface length 

L ss =270.32 mm 

Cascade angle 

y=55 deg 

Cascade flow coefficient 

0=0.80 

Cascade spacing 

S H = 163 mm 

Inlet air angle to the cascade 

=0 deg 

Exit air angle from the cascade 

a 2 = 90 deg 

Rod diameter 

D r = 2.0 mm 

Rod distance to lead, edge 

L r = 122 mm 

Cluster 1 (no rod, steady) 

Sr = 4 00 

O — parameter steady case 

H=0.0 

Cluster 2 rod spacing 

S K = 160.0 mm 

fl — parameter for cluster 1 

fl= 1.59 

Cluster 3 rod spacing 

5^ =80.0 mm 

fl — parameter for cluster 2 

H=3.18 


A computer-controlled three-axes traversing system (Fig. 3) is 
used to measure the inlet velocities, turbulence intensity, and the 
entire boundary layer distribution on suction and pressure sur- 
faces. The traversing system was added very recently to the test 
section to allow the probe to reach all streamwise positions along 
the suction and pressure surfaces, with a positioning accuracy of 
1 fxm. The three-axis system is vertically mounted on the plexi- 
glass side wall. Each axis is connected to a direct current-stepper 
motor with an encoder and decoder. The optical encoder provides 
a continuous feedback to the stepper motor for accurate position- 
ing of the probes. The x- and y-axis of the system are capable of 
traversing along the suction and pressure surfaces in small steps 
up to 1 juum. The third axis rotates the probe holder with an angu- 
lar accuracy less then A® =0.05 deg, which is specifically re- 
quired for boundary layer investigations where the measurement 
of the laminar sublayer is of particular interest. 



Fig. 3 Turbine cascade research facility with three-axis tra- 
versing system 
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Instrumentation, Data Acquisition, and Data Reduction 

The data acquisition system is controlled by a personal com- 
puter that includes a 16-channel, 12-bit analog-digital (A/D) 
board. Time-dependent velocity signals are obtained by using a 
commercial three-channel, constant temperature hot wire an- 
emometer system that has a signal conditioner with a variable low 
pass filter and adjustable gain. A hot wire probe placed upstream 
of the diffuser monitors the reference velocity at a fixed location. 
The pneumatic probes are connected to high precision differential 
pressure transducers for digital readout. Several calibrated ther- 
mocouples are placed downstream of the test section to constantly 
monitor the flow temperature. The wake generator speed and the 
passing frequency signals of the rods are transmitted by a fiber- 
optic trigger sensor. The passage signals of the rods are detected 
by the sensor using a silver-coated reflective paint on one of the 
belts. This sensor gives an accurate readout of the speed of the 
wake generator and the passing frequency of the rods. The signals 
of the pressure transducers, thermocouples and trigger sensors are 
transmitted to the A/D board and are sampled by the computer. To 
ensure the cascade periodicity, the second and fourth blades are 
each instrumented with 48 static pressure taps. Two adjacent 
blades are used for boundary layer measurement. The taps are 
connected to a scani valve, which sequentially transfers the pres- 
sure signals to one of the transducers connected to the A/D board. 

The unsteady data are taken by calibrated, custom designed 
miniature single hot wire probes. At each boundary layer position, 
samples were taken at a rate of 20 kHz for each 100 revolutions of 
the wake generator. The data were ensemble averaged with re- 
spect to the rotational period of the wake generator. Before final 
data were taken, the number of samples per revolution and the 
total number of revolutions were varied to determine the optimum 
settings for convergence of the ensemble average. 

For the steady state case, the instantaneous velocity compo- 
nents are calculated from the temperature-compensated instanta- 
neous voltages by using the calibration coefficients. The instanta- 
neous velocity can be represented in the following form: 

V=V+v (1) 

where V is the mean (time-averaged) velocity and v is the turbu- 
lent fluctuation component. The mean velocity, also known as the 
time average, is given by 



where M is the total number of samples at one boundary layer 
location. The root mean square value of the turbulent velocity 
fluctuation is 
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Re=1 10,000, Tu=1. 9% 


V M 

^g(V y -V) 2 0) 

and the local turbulence intensity is defined as: 


Tuiuc = X 100 = ( Vj - V) 2 X 100 (4) 

V V V M j=l 

For unsteady cases, the ensemble- averaged velocity, fluctuation 
velocity, and the turbulence intensity were calculated from the 
instantaneous velocity samples by: 


= (5) 

7=1 



where A=100 is the total number of wake generator periods and 
M is the number of samples taken per period. (V,(/*)) is the 
ensemble-averaged velocity for the particular boundary layer 
traverse. 

Experimental Results and Discussion 

Detailed surface pressure and boundary layer measurements 
were performed at Reynolds numbers of 110,000 and 150,000. 
These Reynolds numbers, which pertain to a typical cruise opera- 
tion, exhibit a representative value within LPT operating range 
between 75,000 and 400,000, as discussed by Hourmouziadis 
[24]. Furthermore, it produces separation bubbles that can be ac- 
curately measured by miniature hot wire probes. For the Reynolds 
numbers of 110,000 and 150,000, three different reduced frequen- 
cies were examined. To generate unsteady wakes, cylindrical rods 
with the diameter d R = 2 mm were chosen to fulfill the similar 
criterion that requires the generation of a drag coefficient, C D , that 
is approximately equal to the C D of the turbine blade with the 
chord and spacing given in Table 1 (for details we refer to the 
studies in Refs. [23,25]). Furthermore, we define a reduced fre- 
quency f2 that includes the cascade solidity cr, the flow coefficient 
(p, the blade spacing S B , and the rod spacing S R . The reduced 
frequency 12 is an extension of Strouhal number, in the sense that 
it incorporates the rod spacing S R and the blade spacing S B in 
addition to the inlet velocity and wake generator speed. For sur- 
face pressure measurement, rods with uniform spacing as speci- 
fied in Table 1 were attached over the entire belt length. For 
boundary layer measurement, however, clusters of rods were at- 
tached, as mentioned previously. 

Surface Pressure Distributions. Detailed pneumatic surface 
pressure measurements were taken at Re =110,000 and 150,000. 
For each Reynolds number, three different reduced frequencies, 
namely 12=0.0, 1.59, and 3.18, are applied that correspond to the 
rod spacings S R = °° (no rod), 160, and 80 mm. The pressure dis- 
tributions in Fig. 4 show the results of the steady case and two 
unsteady cases. The pressure signals inherently signify the time- 
averaged pressure, because of the internal pneumatic damping ef- 
fect of the connecting pipes to the transducer. The noticeable de- 
viation in pressure distribution between the steady and unsteady 
cases, especially on the suction surface, is due to the drag forces 
caused by the moving rods. The drag forces are imposed on the 
main stream and cause momentum deficiency that leads to a re- 
duction of the total and static pressure. The time-averaged pres- 
sure coefficients along the pressure and suction surfaces are plot- 
ted in Fig. 4. The suction surface (upper portion) exhibits a strong 
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Re=1 50,000, Tu=1 .9% 
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Fig. 4 Static pressure distributions at two different Re num- 
bers and reduced frequencies 12=0, 1.59, 3.18 (no rod, 160, 
80 mm), SS= separation start, SE= separation end 


negative pressure gradient. The flow accelerates at a relatively 
steep rate and reaches its maximum surface velocity that corre- 
sponds to the minimum C p =- 3.5 at s/s 0 =0.42. Passing through 
the minimum pressure, the fluid particles within the boundary 
layer encounter a positive pressure gradient that causes a sharp 
deceleration up to s/s 0 =0.55. This point signifies the beginning of 
the laminar boundary layer separation and the onset of a separa- 
tion bubble. As seen in the subsequent boundary layer discussion, 
the separation bubble characterized by a constant C p plateau ex- 
tends up to s/s 0 ~0.74, thus occupying more than 19% of the 
suction surface and constituting a large separation. Passing the 
plateau, the flow first experiences a second sharp deceleration 
indicative of a process of reattachment, followed by a further 
deceleration at a moderate rate. On the pressure surface, the flow 
accelerates at a very slow rate, reaches a minimum pressure coef- 
ficient at 5/^0 =0.42, and continues to accelerate until the trailing 
edge has been reached. Unlike the suction surface, the pressure 
surface boundary layer does not encounter any adverse positive 
pressure gradient that triggers separation. However, close to the 
leading edge, a small plateau extending from s/s 0 =0.08-0.16 in- 
dicates the existence of a small size separation bubble that might 
be attributed to a minor inlet flow incidence angle. 
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Considering the unsteady case with the reduced frequency XI 
= 1.59 corresponding to a rod spacing of 5^=160 mm, Fig. 4 ex- 
hibits a slight difference in the pressure distribution between the 
steady and unsteady cases. As mentioned above, this deviation is 
attributed to the momentum deficiency that leads to a reduction of 
the total and static pressure. For Re =110,000, the wakes have a 
reducing impact on the streamwise extent of the separation pla- 
teau. As seen in Fig. 4(a), the trailing edge of the plateau has 
shifted from s/s 0 =0.74 to s/s 0 =0.702. This shift reduced the 
streamwise extent of the separation plateau from 19% to 15% of 
the suction surface length which is, in this particular case, 21% of 
reduction in streamwise extent of the separation. Increasing the 
reduced frequency to 11=3.18 by reducing the rod spacing to S R 
= 80 mm causes a slight shift of the C p distribution compared with 
XI =1.59 case. One should bear in mind that pneumatically mea- 
sured surface pressure distribution represents a time integral of the 
pressure events only. 

Increasing the Reynolds number to Re =150,000 has not 
brought major changes in steady state C p distribution. However, 
the combination of higher Reynolds number with unsteady wakes 
reveals the noticeable deviation on the streamwise extent of the 
separation plateau. As seen in Fig. 4(b ) , the trailing edge of the 
plateau has shifted from s/s 0 =0.14 to s/s 0 = 0.702 for Reynolds 
number of 150,000. The combination of higher Reynolds number 
with high unsteady wakes introduce fluctuation kinetic energy into 
the boundary layer which tends to inhibit the separation tendency. 
C p distribution clearly shows that the wake impingement with 
higher Reynolds number shortens the streamwise extent of the 
separation zone compared to the steady case. 

Detailed information regarding the structure of the separation 
bubble is delivered by means of a detailed unsteady boundary 
layer measurement using hot wire probes, as will be discussed in 
the subsequent sections. 

Time Averaged Velocity and Fluctuation Distributions. Fol- 
lowing the surface pressure investigations that mainly addressed 


the onset and extent of the separation zone discussed previously, 
comprehensive boundary layer measurements were performed to 
identify the streamwise and normal extent, as well as the defor- 
mation of the separation zone under unsteady wake flow. The 
steady state case serves as the reference configuration. 

Consistent with the surface pressure distribution above, the ef- 
fect of the wake frequency on the time- averaged velocity profiles 
and fluctuation velocity distribution are presented for one steady 
and two unsteady inlet flow conditions on the suction surface 
along 31 streamwise locations for the Reynolds number of 

110.000, and 41 streamwise locations for the Reynolds number of 

150.000. After completing the velocity measurements, the bound- 
ary layer coordinates were transformed into a blade orthogonal 
coordinate system. Velocities at blade normal positions were ob- 
tained by interpolating their transformed values. The results 
showed almost no difference between the interpolated and nonin- 
terpolated velocity data. Experimental investigations were per- 
formed for three (Afferent values of X1=0.0, 1.59, and 3.18. These 
values cover the reduced frequency range encountered in LPT 
design and off-design operation conditions. 

The effect of wake frequency on time- averaged velocity and 
velocity fluctuation distributions is shown in Figs. 5-8 at six rep- 
resentative streamwise locations for Re =110,000 and Re 
= 150,000. Among 31 streamwise positions measured along the 
suction surface, Figs. 5 and 6 display the velocity and fluctuation 
distributions at one streamwise position upstream, three positions 
within and two positions downstream of the separation bubble. 
The diagrams include the steady state data for reference purposes, 
XI =0.0 (S R =° o), unsteady data for Xl= 1.59 (S R =) 60 mm) and 
Xl=3.18 (5^=80 mm). 

As Fig. 5 indicates, in the upstream region of the separation 
bubble the flow is fully attached. At s/s 0 =0.49, the velocity dis- 
tributions inside and outside the boundary layer experience a 
hardly noticeable decrease with increasing the reduced frequency. 
At the same positions, however, the time-averaged velocity fluc- 
tuations shown in Fig. 6 exhibit substantial changes within the 



Re=110,000, Free stream turbulence with 1.9% 
( without grid), s/s o =0.49 


Re=110,000, Free stream turbulence with 1.9% 
( without grid). s/s„=0.57 


Re=1 10,000, Free stream turbulence with 1 .9% 
( without grid). s/s o =0.61 





Re=1 1 0,000, Free stream turbulence with 1 .9% Re=1 1 0,000, Free stream turbulence with 1 .9% Re=1 1 0,000, Free stream turbulence with 1.9% 
( without grid), s/s o =0.73 ( without grid), s/s 0 =0.77 ( without grid), s/s o =0 85 




Fig. 5 Distribution of time-averaged velocities along the suc- 
tion surface for steady case X1=0 (S R = <») and unsteady cases 
Xl=1.59 (Sfl=160mm) and Xl=3.18 (S fl =80 mm) at Re 
= 110,000 


Re=1 10,000, Free stream turbulence with 1.9% Re=110,000, Free stream turbulence with 1.9% Re=110,000, Free stream turbulence with 1.9% 

( without grid), s/s o =0.73 ( without grid). s/s o =0.77 ( without grid), s/s o =0.85 



Fig. 6 Distribution of time-averaged velocity fluctuations 
along the suction surface for steady case X1=0 (S R =o c) and 
unsteady cases Xl=1.59 (S A =160 mm) and Xl=3.18 (S R 
=80 mm) at Re= 11 0,000 
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boundary layer, as well as outside it. The introduction of the pe- 
riodic unsteady wakes with highly turbulent vortical cores and the 
subsequent mixing has systematically increased the free-stream 
turbulence level from 1.9% in steady case to almost 3% for ft 
=3.18 (5^ = 80 mm). Comparing the unsteady cases, ft =1.59 and 
3.18, with the steady reference case, ft =0.0, indicates that with 
increasing ft the lateral position of the maximum velocity fluc- 
tuation shifts away from the wall. This is due to the periodic 
disturbance of the stable laminar boundary layer upstream of the 
separation bubble. 

As Fig. 6 shows, a substantial influence of the wake frequency 
is observed inside the separation bubble at sls 0 =0.51 , sls 0 
=0.61, and sls 0 =0.13. The wake impingement introduces fluctua- 
tion kinetic energy stemming from its vortical core into the 
boundary layer, trying to energize it and to reverse the separation 
tendency. As seen from the velocity profiles, the wake frequency 
effect shortens the bubble height and reduces its streamwise ex- 
tent. Compared to the steady case, however, the onset of the sepa- 
ration bubble has not changed substantially. This shows that, al- 
though the impingement of the vortical wake core periodically 
reduces the separation bubble height, it does not have sufficient 
momentum to completely suppress it. Figure 6 displays the details 
of the turbulence fluctuation activities in and outside of the sepa- 
ration bubble at s/s o =0.51, s/s 0 = 0.61, and s!s o =0.13. Strong 
turbulent activities are measured within the separation bubble that 
extends to the shear layer. Outside the separation bubble, the tur- 
bulence activities rapidly decrease, approaching the free-stream 
levels that correspond to the individual unsteady frequency dis- 
cussed above. A comparison of the unsteady cases, ft =1.59 and 
3.18, with the steady reference case, ft =0.0, indicates that with 
increasing ft the lateral position of the maximum fluctuation shifts 
toward the wall. This is the consequence of the overall decrease of 
the lateral extension of the separation bubble as a result of wake 
impingement. 

Close to the bubble trailing edge at sls 0 = 0.73 and si s 0 =Q.ll y 
Fig. 5 shows that the undisturbed flow, ft=0 (5^=°o), still has its 


inflectional pattern. The perturbation by impinging the wakes 
causes the velocity profile to become fully turbulent and attached. 
Once the profile becomes fully turbulent, increasing the reduced 
frequency from ft =1.59 (5^ = 160 mm) to ft =3.1 8 (5^ = 80 mm) 
does not substantially change the velocity distribution pattern. 
Figure 6 displays the details of the velocity fluctuation activities 
inside and outside of the reattached boundary layer at si s o =0.11 
and j/5 o = 0.85. As in Fig. 5, once the velocity profile has become 
fully turbulent, its corresponding fluctuation distribution remains 
almost unchanged. 

In the context of transitional and turbulent boundary layer flow 
investigations, it is worth noting the different behavior of these 
boundary layer types. According to the previous investigations by 
Schobeiri et al. [10-12] in a HP- turbine cascade with transitional 
boundary layer throughout, an increased wake frequency causes 
turbulence fluctuations to rise inside and outside the boundary 
layer. However, in the LPT case with the boundary layer separa- 
tion, once the boundary layer is reattached and the velocity distri- 
bution assumes a fully turbulent profile, no major changes are 
observed either in the velocity or in the fluctuation velocity 
distribution. 

For Re =150,000, Figs. 7 and 8 display the velocity and veloc- 
ity fluctuation distributions at the same streamwise locations as 
discussed above. Likewise, the diagrams include the steady state, 
ft =0.0, data for reference purposes and the unsteady data for ft 
= 1.59 and ft=3.18. Although the results for Re= 150,000 show 
very similar velocity and fluctuation patterns, it is observed that 
the starting point of the separation bubble and its reattachment 
point have moved slightly further downstream. This is clearly vis- 
ible in Fig. 7, at s!s o =0.13 and sIs o =0.11 y where the inflectional 
pattern of the velocity profiles indicates that the above streamwise 
positions are located close to the trailing edge within the separa- 
tion bubble. Increasing the frequency from ft=0 (S R = °°) to ft 
=3.18 (5^=80 mm) brings the profiles closer to the fully turbulent 
pattern, without collapsing. The same statement is true for the 
turbulence fluctuations, Fig. 8, at sls a = 0.73 and s!s 0 =Q.ll. Also, 


Re= 150,000, Free stream turbulence with 1.9% 
( without grid), s/s o =0.49 



Re=1 50,000, Free stream turbulence with 1.9% 
( without grid), s/s ,.0,73 



Re=1 50,000, Free stream turbulence with 1.9% Re= 150,000, Free stream turbulence with 1.9% 
( without grid). s/s o =0.57 ( without grid), s/s o =0.61 




Re=1 50,000, Free stream turbulence with 1 .9% Re=1 50,000, Free stream turbulence with 1.9% 

( without grid), s/s o =0.77 ( without grid), s/s„=0.85 




Re= 150.000, Free stream turbulence with 1.9% 
( without grid), s/s o =0.49 



Re=1 50,000, Free stream turbulence with 1 .9% Re=1 50,000, Free stream turbulence with 1.9% 

( without grid), s/s o =0.57 ( without grid), s/s o =0.61 



Re=150,000, Free stream turbulence with 1.9% Re=1 50,000, Free stream turbulence with 1.9% Re=150,000, Free stream turbulence with 1.9% 

( without grid), s/s o =0. 73 ( without grid), s/s„=0.77 ( without grid), s/s o =0.85 



Fig. 7 Distribution of time-averaged velocity distributions 
along the suction surface for steady case ft=0 (S R = «) and 
unsteady cases 11=1.59 (S fl =160 mm) and 11=3.18 (S fl 
=80 mm) at Re =150,000 


Fig. 8 Distribution of time-averaged fluctuation velocities 
along the suction surface for steady case 11=0 (S R = °o) and 
unsteady cases 11=1.59 (S fl =160 mm) and 11=3.18 ( S R 
=80 mm) at Re =150,000 
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the size of the separation bubble is smaller when compared to that 
for Re= 11 0,000. To achieve a noticeable impact of the Reynolds 
number on the onset and extent of the bubble, the Reynolds num- 
ber needed to increase above 150,000. 

Temporal Behavior of the Separation Zone Under Periodic 
Unsteady Wake Flow. Velocity distributions on the suction sur- 
face with time as the parameter are plotted in Figs. 9-12 for Re 
= 1 10,000, 150,000 for rod spacings S*=80 mm and S R 
= 160 mm that correspond to 11=1.59 and 3.18. The nondimen- 
sional time (t/ r) values are chosen that they represent the tempo- 
ral states within one full period of wake passing. For Re 
= 110,000, Figs. 9(a)-9(f) show the velocity distributions, inside 
and outside the boundary layer at fixed s/s Q locations, experience 
moderate to pronounced changes. Figure 9(a) represents the in- 
stantaneous velocity distribution upstream of the separation zone, 
followed by Figs. 9(b)-9(f) which represent the velocity distribu- 
tions inside the separation zone. In discussing the following re- 
sults, we simultaneously refer to the wake distribution as well as 
the turbulence fluctuation results. 

Figure 9(a) exhibits the velocity distribution on the suction sur- 
face at s!s o =0A9. At this streamwise position, the laminar bound- 
ary layer is subjected to a strong negative pressure gradient. The 
velocity distributions at different (t! t) experience changes in 
magnitude that reflect the corresponding changes of the impinging 
periodic wake velocity. It is worth noting that, despite the injec- 
tion of turbulence kinetic energy by the impinging wakes, no local 
instantaneous boundary layer transition occurs. This is because of 
the strong negative pressure gradient that prevents the boundary 
layer from becoming instantaneously transitional. Instantaneous 
velocity distributions inside the separation zone are shown in Figs. 
9(b)-9(f). 

As a representative case, we discuss the results plotted in Fig. 
9(e) at s/s o =0.65. During the time interval from tlr close to 0.5 
(1.5, 2.5, etc.) to about tlr-0.15 (1.75, 2.75, etc.), the separation 
zone is exposed to the wake external flow with relatively lower 


turbulence level. This flow does not have the capability to sup- 
press the separation zone. Thus, the separation region is clearly 
shown by the velocity distributions at tlr= 0.5 and //t= 0.75. As 
the wake passes over the blade at sls Q = 0.65 introducing high 
turbulence kinetic energy into the boundary layer, the boundary 
layer is energized, causing the separation zone to partially reduce. 
To emphasize this statement, the steady state velocity distribution 
at the same streamwise position is also plotted in Fig. 9(e) using 
full circles. It shows clearly the separated nature of the boundary 
layer which coincides with the instantaneous velocity profile at 
tlr=QJ5. Intermediate times reflect the gradual change between 
different tlr states as the flow is undergoing the influence of the 
oncoming wake. Moving to the trailing edge of the separation 
zone, at s/s^O.73, Fig. 9(f ) , a partial reduction in boundary layer 
thickness as the result of wake impingement is visible, however, 
the separation zone does not seem to disappear. 

Figures 10(a)-10(/) show that the velocity distributions outside 
the boundary layer at fixed sls 0 locations experience noticeable 
changes at 11=3.18 (5^ = 80 mm). Increasing the wake passing 
frequency causes the wake turbulence kinetic energy to increase, 
resulting in a stronger suppression compared with the 11 = 1.59 
(5^=160 mm) case. As seen in Figs. 11 and 12, increasing the 
Reynolds number to 150,000 causes the leading and trailing edges 
of the separation bubble to move further downstream. It is also 
observed that increasing the Reynolds number reduces the size of 
the separation bubble. 

The effect of the periodic unsteady wakes on the onset and 
extent of the separation bubble is shown in Figs. 13-16 for Re 
= 110,000 and 150,000, and for two different frequencies, namely 
0=1.59 and 0=3.18. These figures display the full extent of the 
separation bubble and its behavior under a periodic wake flow 
impingement at different tlr. The wake propagation for 0 = 1.59 
and 0=3.18 is analyzed and the value of t! r corresponds to the 
point in the cycle at which the data acquisition system is triggered. 
During a rod passing period, the wake flow and the separation 
bubble undergo a sequence of flow states which are not noticeably 


Re= 11 0,000, Free stream turbulence with 1.9% 
( without grid), Q=1 .59, s/s o =0.49 


Re=1 1 0,000, Free stream turbulence with 1 .9% 
( without grid), n=1.59, s/s o =0.53 


Re=1 1 0,000, Free stream turbulence with 1 9% 
( without grid), 0=1 .59, s/s o =0.57 


Re= 11 0,000, Free stream turbulence with 1.9% 
( without grid), n=3.1 8, s/s o =0.49 


Re=1 10,000, Free stream turbulence with 1 .9% 
( without grid). 0=3.1 8. s/s 0 =0.53 


Re=1 10,000, Free stream turbulence with 1.9% 
( without grid), 0=3.18, s/s o =0.57 






V(m/s) V(m/s) V(m/s) 


Re=1 1 0,000, Free stream turbulence with 1 .9% 
( without grid), 0=1 .59, s/s„=0.61 


Re=1 1 0,000, Free stream turbulence with 1 .9% 
( without grid), 0=1.59, s/s o =0.65 


Re=1 1 0,000, Free stream turbulence with 1 .9% 
( without grid), 0=1 .59, s/s o =0 73 


Re=110,000, Free stream turbulence with 1.9% 
(without grid), 0=3.18, s/s o =0.61 


Re=1 1 0,000, Free stream turbulence with 1 .9% 
( without grid). 0=3. 1 8. s/s o =0.65 


Re= 11 0,000, Free stream turbulence with 1.9% 
( without grid), 0=3.18, s/s„=0.73 
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Fig. 9 Distribution of the ensemble-averaged velocity devel- Fig. 10 Distribution of the ensemble-averaged velocity devel- 
opment along the suction surface for different s/Sq with time opment along the suction surface for different s/Sq with time 

f/r as parameter for 11=1.59 (S fl =160 mm) and Re=1 10,000 f/r as parameter for 11=3.18 (S R = 80 mm) and Re=1 10,000 
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Re=150,000, Free stream turbulence with 1.9% Re=1 50,000, Free stream turbulence with 1.9% Re=150,000, Free stream turbulence with 1.9% 

( without grid), fi=1 .59, s/s o =0.49 ( without grid), 0=1 .59, s/s o =0.53 ( without grid), 0=1 .59, s/s o =0.57 



V(m/s) V(m/s) V(m/s) 


Re=150,000, Free stream turbulence with 1.9% Re=150,000, Free stream turbulence with 1,9% Re=150,000, Free stream turbulence with 1.9% 

( without grid), 0=1 .59, s/s o =0.61 ( without grid), 0=1 .59, s/s o =0.65 ( without grid), 0=1 .59, s/s„=0.73 



V(m/s) V(m/s) V(m/s) 

Fig. 11 Distribution of the ensemble-averaged velocity devel- 
opment along the suction surface for different s/Sq with time 
f/ras parameter for fl=1.59 (S fl =160 mm) and Re=150,000 



V(m/s) V(m/s) V(m/s) 


Re=1 50,000, Free stream turbulence with 1 .9% Re=1 50,000, Free stream turbulence with 1 .9% Re=1 50,000, Free stream turbulence with 1 .9% 

( without grid), 0=3.18, s/s„=0.61 ( without grid), 0=3.18, s/s„=0.65 ( without grid), 0=3.18, s/s o =0.73 



V(m/s) V(m/s) V(m/s) 


Fig. 12 Distribution of the ensemble-averaged velocity devel- 
opment along the suction surface for different s/Sq with time 
tlr as parameter for ft =3.1 8 (S fl =80 mm) and Re= 150,000 


different when the unsteady data are time averaged. 

Starting with Re =100,000 and ft =1.59, Fig. 13(a) exhibits the 
separation bubble in its full size at t/r= 0.25. At this instant of 
time, the incoming wakes have not reached the separation bubble. 
At tlr= 0.5, the wake, with its highly turbulent vortical core, 
passes over the blade and generates high turbulence kinetic en- 
ergy. At this point, the wake turbulence penetrates into the bubble 
causing a strong mass, momentum, and energy exchange between 
the wake flow and the fluid contained within the bubble. This 
exchange causes a dynamic suppression and a subsequent contrac- 
tion of the bubble. As the wake travels over the bubble, the size of 
the bubble continues to contract at t/r=0J5 and reaches its mini- 
mum size at t / r= 1 .0. At // r= 1 , the full effect of the wake on the 
boundary layer can be seen before another wake appears and the 
bubble moves back to the original position. Similar results are 
observed when operating at the same Reynolds number, Re 
= 110,000, but at a higher reduced frequency, ft =3. 18, Figs. 
14(a)-14(d). A comparison of Figs. 14(a)-14(d) with Figs. 
\2{a)-\2{d) illustrates the pronounced influence of higher reduced 
frequency which causes more frequent penetration of wake turbu- 
lence into the separation bubble, causing a stronger suppression 
and subsequent contraction of the bubble. 

Figures 15 and 16 display the behavior of the separation bubble 
under higher Reynolds number, Re =150,000, and reduced fre- 
quencies, ft =1.59 and ft =3. 18. While the higher Reynolds num- 
ber moves the separation bubble further downstream. Similar con- 
tracting effects as discussed above are observed at ft =1.59. 
However, doubling the reduced frequency and increasing the Rey- 
nolds number is associated with the higher turbulence intensity 
that leads to stronger suppression of the separation bubble, as 
shown in Figs. 15 and 16. 

Intermittency Analysis 

Intermittency distribution, which identifies whether the flow is 
laminar or turbulent inside the boundary layer, is calculated fol- 
lowing the method of Hedley and Keffer [26]. Instantaneous ve- 


locity is sensitized to increase its discriminatory capabilities be- 
tween turbulent and nonturbulent parts of the signal. For this 
purpose, the multiplication of the first derivative of the velocity 
signal and the velocity signal is used for further analysis. This is 
called the detector function, 5(t), defined as 


S(t) = 



(8) 


The above detector function is used by many researchers, includ- 
ing Antonia and Bradshaw [27], Kovasznay et al. [28], Bradshaw 
and Murlis [29], and Schobeiri and his coworkers [7,12]. A very 
recent paper by Schobeiri [30] extensively discusses different is- 
sues relative to intermittency based transition modeling. Though 
sensitized detector function separates the turbulent and nonturbu- 
lent zones of the fluid, there is still some overlap between the two 
near the origins. The discrimination between the two zones of the 
flow will be ideal when the overlap between the two distributions 
is minimal or zero. To eliminate the disturbing effects of the ve- 
locity signal peaks, a smoothing procedure is applied to the S(t) 
signal. The mean value of ten consecutive S(t) values is calculated 
and the ten values are substituted by their mean value of S sm (t). 
After smoothing the detector function, a threshold level C is then 
applied to the smoothed detector function to distinguish between 
true turbulence and the signal noise. It is defined as 


1 1 when S sm (t) 5* C 
[o when S sm (t) < C 


(9) 


Once the threshold level is applied to the detector function S(t), 
the result is a random square wave called the indicator function /, 
with 0’s representing the nonturbulent case and l’s indicating the 
turbulent behavior of the boundary layer. A threshold level, C, of 
1.2 is used for all the data on the suction surface. In the absence of 
length scales, this value is chosen from visual observations. Sev- 
eral other values of C are tested and little qualitative difference is 
seen in the intermittency distribution during transition. Though the 
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s/s D s/s Q 

Fig. 13 Ensemble-averaged velocity contours along the suction surface for different s/Sq with 
time tlr as parameter for 11=1.59 (S fl =160 mm), Re=110,000 


(a) Re-1 10,000, 0-3.18, Vx-0.25 (b) Re-110,000, Q-3.18, t/x-0.5 



s/s Q s/s D 

Fig. 14 Ensemble-averaged velocity contours along the suction surface for different s/Sq with 
time f/r as parameter for n=3.18 (S fl =80 mm), Re= 11 0,000 


intermittency values vary with different values of C, the important 
parameters like start and end of transition are not affected by C. 
The resulting square wave after applying the threshold is en- 
semble averaged to get the ensemble-averaged intermittency as 
follows 

Journal of Turbomachinery 


= ( 10 ) 

™j = 1 

where N is the number of revolutions of the wake generator for 
which the data are collected. For time- averaged intermittency, 
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Fig. 15 Ensemble-averaged velocity contours along the suction surface for different s/Sq with 
time tlr as parameter for 11=1.59 (S fl =160 mm), Re=150,000 


(a) Re-150,000, 0=3.18, t/x-0.25 (b) Re-1 50,000, 0-3.18, t/i-0.50 



Fig. 1 6 Ensemble-averaged velocity contours along the suction surface for different s/Sq with 
time tlr as parameter for 11=3.18 {S R =S0 mm), Re= 150,000 


(y f (^)) is integrated with respect to time to arrive at 



Figure 17 shows the processing of instantaneous velocities. 

102 / Vol. 129, JANUARY 2007 


Ensemble-Averaged Intermittency Distribution. Figure 18 
presents the temporal-spatial contours of the ensemble-averaged 
intermittency distribution at three different lateral positions above 
the blade suction surface, y= 1.341, 1.755, and 6.10 mm, for two 
reduced frequencies, 11=1.59 and 11=3.18. For better comparison 
of the effects of the impinging wake frequency, Fig. 18 exhibits 
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only the first three wakes. In this figure, the wakes with the highly 
vortical cores display intermittency values close to 0.6 indicating 
the transitional character of the boundary layer at the particular 
instant of time of wake impingement on the surface. Intermittency 
is approximately equal to zero outside the wake region near the 
leading edge, showing the nonturbulent behavior of the flow. 

The wakes represented by narrow green strips pass through the 
turbine blade channel and periodically switch the boundary layer 
from laminar to turbulent and vice versa. Upstream of the separa- 
tion bubble a pretransitional strip with an ensemble-averaged in- 
termittency of y(/)~0.5 stalling at £/£ 0 ~0.43 and ending at 
s/ S q ~ 0.52 separates the attached boundary layer from the bubble 
leading edge. At s/s 0 ~0.52, the visibility of the wake vanishes 
due to the interaction with the separation bubble. As Fig. 18 
shows, the separation bubble starts at s/s 0 ~0.52 and extends up 
to s/sq~ 0.75, thus occupying more than 24% of the suction sur- 
face and forming a large separation zone. At s/s 0 ~OJ5 , the in- 
termittency field in Figure 18(a) displays an abrupt change in 
intermittency level, which indicates the start of a reattachment 
process. Once the wake passes over the separation zone, its sig- 
nature reappears again as spots, with higher intermittency level 
(red) associated with becalmed zones (blue). Increasing II to 3.18, 
Fig. 18(6), causes an earlier mixing of the impinging wakes which 
results in widening the areas occupied by the wake vortical core, 
thus reducing the wake external region. Figure 18(6) suggests that 
further increase of 11 may lead to a complete degeneration of the 
deterministic periodic wake flow into a stochastic turbulence. Fig- 
ures 18(c)-18(/) display the intermittency distributions at higher 
normal position from the blade surface. They are quite identical 
with Figs. 18(a) and 18(6) discussed above. 

The intermittency distributions in Fig. 18 clearly show the un- 
steady nature of the boundary layer transition. In this form, how- 
ever, they cannot quantitatively describe the complex unsteady 
transition, separation and reattachment process. To establish the 
basic relations essential for a quantitative description of the un- 
steady boundary layer transition, we resort to the fundamental 
study by Schobeiri and his coworkers [23] that deals with the 
physics of steady and unsteady wake development in a curved 
environment. The study clearly shows that the turbulence structure 
of the steady and unsteady wake flow is determined by the wake 
defect, which is a Gaussian function. Following the above study, 
we define a dimensionless parameter 

/»+oo 

c = =t ^ = ~h with b= ~r\ rd & ( 12 ) 

b ib b Vtj- J_ oc 

Equation (12) relates the wake passing time t with the wake pass- 
ing velocity in the lateral direction U w , and the intermittency 
width 6, with as the lateral distance from the wake center [23]. 
The intermittency width 6 is directly related to the wake width 
introduced by Schobeiri and his coworkers [7,12,23]. In an analo- 
gous way to find the defect function, we define the relative inter- 
mittency, T, as 


"(r^)L-(rA)L 1 

In the above equation, (%•(**)) the time dependent ensemble- 
averaged intermittency function, which determines the transitional 
nature of an unsteady boundary layer. The intermittency (y,(^))max 
exhibits the maximum intermittency value inside the wake vorti- 
cal core. Finally, (y,(^))min represents the minimum ensemble- 
averaged intermittency values outside the wake vortical core. 

A representative relative intermittency function, T, is shown in 
Figs. 19(a)-19(d) for a frequency value of 11=1.59 at lateral dis- 
tances from the blade surface of y =0.858, 0.996, 5.3, and 9.3 mm, 
with the dimensionless longitudinal distance s/s 0 as a parameter. 
The above distances are representative for intermittency distribu- 
tions before, inside and outside the separation bubble over the 


entire suction surface. The symbols represent the experimental 
data. As seen for the reduced frequency of U= 1.59, the measured 
relative intermittency functions follow very closely a Gaussian 
distribution, given by 


r = e-? (14) 

with f as the nondimensionalized lateral length scale defined in 
Eq. (12). The slight deviation within and after the separation zone 
is due to a difficulty associated with capturing the maximum and 
minimum intermittencies. Using Eq. (14) as a generally valid in- 
termittency relationship for unsteady wake flows, the intermit- 
tency function (y*(fj)) is completely determined if additional in- 
formation about the intermittency functions (%•(*;)) max and 
(y*W)min are available. The distribution of (y f (f f )) max and 
(%•(/*)) min i n the streamwise direction are plotted in Fig. 20 for II 
values of 1.59 and 3.18 on the suction surface. The distribution of 
( %•(£{) max) corresponds to the condition when the wake, with its 
high turbulence intensity core, impinges on the plate surface. 
Once the wake has passed over the surface, the same streamwise 
location is exposed to a low turbulence intensity flow regime with 
an intermittency state of (‘y l (f f )) min , where no wake is present. 

Figure 20 displays the striking features of (y(O)max and 
While for zero and moderate pressure gradients, the 
minimum intermittency (%•(*;)) min distribution reveals a certain 



Fig. 17 Calculation of ensemble-averaged intermittency func- 
tion from instantaneous velocities for 11 = 1.725 at y 
=0.720 mm 


Journal of Turbomachinery 


JANUARY 2007, Vol. 129 / 103 


Downloaded 10 Dec 2008 to 128.156.10.80. Redistribution subject to ASME license or copyright; see http://www.asme.org/terms/Terms_Use.cfm 





Fig. 18 Ensemble-averaged intermittency factor in the temporal-spatial domain at different y 
positions for 11=1.59 (S fl =160 mm) and 0=3.18 (S fl =80 mm) 


similarity to the one described by the Emmons-Narasimha transi- 
tion model [31], the present LPT- flow case with a strong negative 
pressure gradient associated with separation, (y,(^))min exhibits a 
remarkably different course that occurs systematically and repro- 
ducibly for all O cases at all y positions over the blade surface. 
Upstream of the separation bubble, the course of (y(O)min with the 
value close to zero indicates a stable nonturbulent character of the 
boundary layer. A sharp increase in intermittency indicates the 
separation begins shortly before the pressure minimum (Fig. 4) 
has been reached. It is followed by a high intermittency region 
that covers the separation plateau (Fig. 4) and a steep decrease 
that is indicative of reattachment. The streamwise location of the 
intermittency minimum at sls 0 ~ 0.7 coincides with the end of the 
separation plateau. The following increase in intermittency is due 
to the steep positive pressure gradient that follows the constant 
pressure plateau. On the other hand, (y,(0)max reveals a funda- 
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mentally different behavior. As Fig. 20 shows, the wake flow with 
an intermittency of (y f (^)) ma x~0. 8-0.9 impinges on the blade 
surface. By convecting downstream, because of an extremely thin 
boundary layer, the wake turbulent fluctuations do not undergo a 
strong damping by the wall shear stress forces, as was observed in 
zero and moderate pressure gradient cases reported in Refs. 
[7,12]. 

Utilizing (y(f))max and (yM)min> the relative intermittency T is 
found to be described by a Gaussian distribution. This observation 
is in accord with the findings reported in Refs. [7,12], and very 
recently in Ref. [30] confirming the universal character of T. Con- 
sidering the intermittency results of the current investigations and 
those reported above, it can be concluded that, in general, 
(y*W)max and (y f -(0)min are not only functions of reduced fre- 
quency, but they are also strongly influenced by the pressure gra- 
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Fig. 19 Relative intermittency as a function of s/Sq for unsteady frequency of 11 = 1.59 ( S R 
= 160 mm) at (a) y= 0.858 mm, (b) y= 0.996 mm, (c) y=5.3 mm, and (d) y=9.3 mm at Re 
= 110,000 


dient, turbulence intensity, Reynolds — and possibly Mach num- 
ber, and surface roughness. This implies that neither (y { (/ { )) m ^ nor 
(%•(*,•)) max have universal character. 

Uncertainty Analysis 

The Kline and McClintock [32] uncertainty analysis method 
was used to determine the uncertainty in the velocity after cali- 
bration and data reduction for the single- wire probe. The Kline 
and McClintock method determines the uncertainty with a 95% 
confidence level. The uncertainty in the velocity for the single- 
wire probe after the data reduction is given in Table 2. As shown, 
the uncertainty in the velocity increases as the flow velocity de- 
creases. This is due to the pneumatic pressure transducer having a 
large uncertainty during calibration. 

Conclusions 

A detailed experimental study on the behavior of the separation 
bubble along the suction surface of a highly loaded LPT blade 
under periodic unsteady wake flow was presented. Varying the 
Reynolds number, one steady and two different unsteady inlet 
wake flow conditions with the corresponding passing frequencies, 
the wake velocity and the turbulence intensities were investigated 
by utilizing a large-scale, subsonic research facility. Periodic un- 
steady wake flow was established by translational motion of two 
parallel moving timing belts on which cylindrical rods were at- 
tached. The following conclusions were drawn. 


Table 2 Uncertainty in velocity measurement for hot-wire 
probe 


V (m/s) 

3 

5 

12 

(O V/V 0 (%) 

5.78 

2.41 

1.40 
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(1) Slight changes of the pressure distribution occurred, while 
operating at the unsteady flow conditions. Increasing the Reynolds 
number from Re= 110,000 to Re =150,000 has not brought major 
changes in steady state C p distribution. However, the combination 
of higher Reynolds number with higher unsteady wake frequency 
introduced higher fluctuation kinetic energy into the boundary 
layer, which tends to reverse the separation tendency. C p distribu- 
tion clearly shows that the wake impingement with higher Rey- 
nolds number shortens the streamwise extent of the separation 
zone compared to the steady case. 

(2) Detailed unsteady boundary layer measurement identified 
the onset and extent of the separation bubble as well as its behav- 
ior under the unsteady wake flow. Passing the wake flow with its 
highly turbulent vortical core over the separation region caused a 
periodic contraction and expansion of the separation bubble and a 
reduction of the separation bubble height. Increasing the passing 
frequency associated with a higher turbulence intensity further 
reduced the separation bubble height. It was observed that, by 
increasing the Reynolds number to 150,000, the leading edge and 
trailing edge of the separation bubble and, thus, the re-attachment 
point, moved further downstream to s/so=0.56 and s/sq =0.788, 
respectively. Also, the size of the separation bubble was further 
reduced. 

(3) Intermittency analysis of the current boundary layer experi- 
mental data with the flow separation determined the minimum, 
maximum, and the relative intermittency functions, (y min ), 
-<y max ), and T. The minimum intermittency function, (y min ), rep- 
resented the boundary layer behavior when it is exposed to the 
wake external region (region between the turbulent wake strips). 
On the other hand, (y max ) describes the state of the boundary layer 
when it is subjected to the wake vortical core with its high turbu- 
lence level. 

(4) The relative intermittency factor followed a Gaussian distri- 
bution confirming the universal character of the relative intermit- 
tency function. In contrast to the relative intermittency function T, 
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0=3.18, y=0.765 mm 



Fig. 20 Maximum, minimum, and time-averaged intermittency as a function of 
s/Sq at different lateral positions for steady case 0=0 (S R = «=) and unsteady cases 
0=1.59 (S R = 160 mm) and 0=3.18 (Sk= 80 mm) at Re= 11 0,000 


the minimum as well as the maximum intermittency do not sug- 
gest to have a universal character. Several parameters, such as 
free-stream turbulence intensity, reduced frequency, surface 
roughness, Re number, and pressure gradient, are instrumental in 
affecting the pattern of these two intermittencies. Future studies 
need to incorporate these parameters. 
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Nomenclature 

b = intermittency wake width 
c = blade chord 
c ax = axial chord 

Cp = pressure coefficient, C p =(pj-p s )/(p t -p s ) M 
C = threshold level 


d R 

H n 

h m 

^ss 
M 
N 
Pi 
Ps > Pt 

Re LSS 


^md 


So 

S r 

s s 

S(x,t) 


rod diameter 

shape factor, H n =S\/ 

maximum separation bubble height 

indicator function 

suction surface length 

number of samples 

number of wake cycles 

static pressure taps i- 1 , . . ,48 

static, total pressure 

suction surface Reynolds number 

— ^ss^exit^ ^ 

streamwise distance from the leading edge of 
the blade 

streamwise location of maximum separation 
bubble height 

streamwise distance from the leading edge to 
the trailing edge of the blade 
re-attachment point of the separation bubble 
from blade leading edge 
start of the separation bubble at a streamwise 
distance from blade leading edge 
criterion function 
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S B = blade spacing 
S R = rod spacing 
t = time 

Tu = reference turbulence intensity 
(Tu) = ensemble- averaged turbulence intensity 
U = belt translational velocity 
v = fluctuation velocity 
V = velocity 
V ax = axial velocity 
V exit — exit velocity 
x = position vector 
y = lateral distance from plate surface 
a 12 — cascade inlet, exit flow angles 
y = cascade stagger angle 
y = time- averaged intermittency 
((y)) = ensemble- averaged intermittency 
<(y)> max maximum ensemble- averaged intermittency 
((y))min = minimum ensemble- averaged intermittency 
T = relative turbulence intermittency 
l — nondimensional coordinate, yib 
v = kinematic viscosity 
£2 = lateral distance from wake center 
p = density of air 
a = cascade solidity, a=c/S B 
7 = one wake-passing period 
(p = flow coefficient, (p=V m IU 
ip A = Zweifel coefficient 

^ 4 = 2 sin 2 a 2 (cot a2~cot a\)S B lc „ 

II = reduced frequency ( 1 =(clS R )(U/ V^) 
=(<r/<p)(S B /S R ) 
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